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Abstract 

We analyse the 30 October, 2004, X1.2/SF solar event that occurred in AR 10691 
(N13 W18) at around 11:44 UT. Observations at 212 and 405 GHz of the Solar 
Submillimeter Telescope (SST), with high time resolution (5 ms), show an intense 
impulsive burst followed by a long-lasting thermal phase. EUV images from the 
Extreme Ultraviolet Imaging Telescope (SOHO/EIT) are used to identify the pos- 
sible emitting sources. Data from the Radio Solar Telescope Network (RSTN) com- 
plement our spectral observations below 15 GHz. During the impulsive phase the 
turnover frequency is above 15.4 GHz. The long-lasting phase is analysed in terms of 
thermal emission and compared with GOES observations. From the ratio between 
the two GOES soft X-ray bands, we derive the temperature and emission measure, 
which is used to estimate the free-free submillimeter flux density. Good temporal 
agreement is found between the estimated and observed profiles, however the former 
is larger than the latter. 
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1 Introduction 



It is widely accepted tliat solar flares imply the release of large amounts of 
energy, which go mainly in the acceleration of particles: electrons and ions. 
The interaction of these particles with the medium, plasma and magnetic fleld, 
produces radiation by very different mechanisms, which is observed at Earth 
in a wide range of wavelengths. Of particular interest are electrons moving 
along the magnetic lines and producing gyrosynchrotron emission observed 
at microwaves and shorter wavelengths. Gyrosynchrotron emission of mildly 
relativistic electrons is distributed in a broad continuum at high harmonics of 
the fundamental frequency ub = eB/{27rmec), where e is the electron charge, 
B the magnetic fleld strength, rrie the electron mass at rest and c the speed of 



light. The maxi mum emissio n is centred around v 
Lorentz factor (Dulk 



1^57 , with 7 the electron 



19851). Taking 100 < B < 1000 G, indi cative values 



of 7 a re between 2 and 5 for a 10 GHz photon. Observationally, iKosugi et al 



( 1l988l ) found that the highest correlation between centimetric wavelengths and 
hard X-rays (HXR) happens between 17 GHz and 80 keV; that is, electrons 
with energies < 200 keV. Therefore, radio observations at shorter wavelengths 
offer new insights to understand processes that involve high energy electrons 
during solar flares, where current HXR detector technology does not provide 
enough sensitivity and/or signal to noise ratios to observe weak flares. Since 
1999, the only submillimeter instrument dedica ted to solar observation is the 
Solar Submillimeter Telescope (SST, see e.g.. 



Kaufmann et al 



20011) oper- 
ating at 212 and 405 GHz. From the theoretical arguments given above, > 4 
MeV electrons should produce the synchrotron emission observed with SST. 



In this work, we report the multiwavelength observations of the impulsive and 
gradual phases of an X-class flare that occurred on 30 October, 2004. In the 
following section we describe the instruments and data analysis. In Section 3, 
we analyse the impulsive phase of the flare; while in Section 4, we do the same 
for the gradual phase. Finally, we give our concluding remarks. 



2 Observations and principal characteristics of the event 



The SST (IKaufmann et al.l . 1200 ll ) is a radome enclosed 1.5 m single dish an- 
tenna with room temperature receivers in a focal array: 4 at 212 GHz and 2 at 
405 GHz, with 5 ms time resolution. On 30 October, 2004, SST was tracking 
AR 10691 (N13 W18) with its beam 5. In Figure 1, we sketch the position of 



SST beams with crosses on a 195 A EUV image of EIT (jPelaboudiniere et al. 



19951 ): the size of the crosses is equal to the HPBW. While beams 2, 3 and 4 



(212 GHz) are round, beam 5 (405 GHz) has an elliptical shape, with minor 
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Fig. 1. A 195 A EIT image of the flaring area at 11:36:42 UT. North is up and 
east to left. The labelled crosses represent the SST beams: 2, 3 and 4 are 212 GHz 
receivers, while 5 is a 405 GHz receiver. Only beams 3 and 5 detected the burst. 



and major axes in the directions shown in the figure. Because of the high at- 
mospheric attenuation only two beams observed this event, one at 212 GHz 
(beam 3) and one at 405 GHz (beam 5). Beams 1 (212 GHz) and 6 (405 GHz) 
are not shown in the figure. Since only two beams ob served the event, we can- 
not use the multibeam technique (j Cost a et al.l . Il995l ) to locate the centroid of 
emission. 



Assuming that the radio source is coincident with the EUV bright loop (see 
Fig. 1), we correct the observed antenna temperature at 212 and 405 GHz 
for the misalignment. During calibration the antenna temperatures are also 
corrected for atmospheric attenuation. Finally, the corrected antenna tem- 
peratures are converted to flux densities, F, using the well known formula: 
F = 2kbTa/A(., where kf, is the Boltzmann constant, the corrected antenna 
temperature and A^, the effective area of the antenna. We estimate a 40% un- 
certainty for the whole process. 



Data at microwave frequencies with 1 second time resolution were obtained 



by the Radio Solar Telescope Network (RSTN. iGuidice et al.l . Il98ll ). Figure 2 



shows the time profiles of the event at different wavelengths. At 15.4 GHz, we 
see that the burst starts at 11:42:40 UT as a small enhancement that evolves 
into the impulsive phase with a maximum flux density around 3300 s.f.u. at 
11:44:09 UT. The increasing phase is characterized by short emissions that can 
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Fig. 2. Time profiles at selected wavelengths. GOES data have 3 seconds time res- 
olution, radio data are integrated to 1 second. 

be attributed to different particle injections. The decreasing phase is smooth, 
with no bursts, reaching the minimum value (close to the pre-fiare level) at 
11:45:25 UT. Afterward, a post burst increase is observed until approximately 
11:51 UT. This gradual phase has also some bursting episodes characteristic 
of particle injections. 

Although the 212 GHz time profile is in general smoother than that at 15.4 
GHz, most of the features observed at 15.4 GHz are also present, with smaller 
amplitude, at 212 GHz. The starting time of the burst at 212 GHz coincides 
with that at 15.4 GHz (11:42:40 UT), but the peak time is delayed by 3 
seconds (11:44:12 UT at 212 GHz). The gradual phase is also present with 
sudden changes, probably, of instrumental origin. The 405 GHz time profile 
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Fig. 3. Wavelet transform of the flux density at 15.4 and 212 GHz and GOES time 
derivative represented in a 2D graph. Contour levels are: 10~^, 10~^ W m~^ s~^ 
(GOES time derivative); 75, 750 s.f.u. (212 GHz) and 185, 1850 s.f.u. (15.4 GHz). 

during the rising phase is closer to the 15.4 GHz emission than to the 212 
GHz, is in general delayed in 6 seconds. The decaying phase is more extended 
than in the other frequencies and shows the presence of longer features, some 
of these features could have an atmospheric origin, and do not allow us to 
clearly recognize the gradual phase. 



Unfortunately, no HXR data are available for this event. Instead, we have 
used the time derivative of GOES 1.5-12 keV soft X-ray (SXR) channel with 
3 seconds temporal resolution to compare with the radio flux density. The 
sensitivity of GOES is better than 5 10~^ W m~^, taking this value as the 
flux uncertainty, the flux time derivative has a mean uncertainty better than 
V2 5 10-73 W m-2 s ~ 3 10-8 W m-^ s-\ The time derivative of GOES 
channel 3 (bottom curve in Fig. 2) mimics very well the radio time profiles 
during the rising phase and is coincident in time, but has two clear maxima, 
the first peaking in coincidence with the 15.4 GHz emission. We will show later 
that the second GOES time derivative peak, still during the rising phase, is 
also present at the other frequencies, but during their decaying phase. Finally, 
the decaying phase in this curve lasts longer than the others and has some 
bursting episodes. 



We use the wavelet transform (see e.g. iDaubechied . Il992l ) to compare the 
time profiles. Wavelets are very useful to discriminate the different time scales 
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in a time series. We have integrated our data to the lowest time resolu- 
tion, namely 3 seconds and applied the wavelet transform implemented by 



Torrence and Compd (119981 ) based on a Morlet mother wavelet with wavenum- 



ber 6. In Figure 3 each contour represents a structure resolved by the transform 
with a characteristic scale in the ordinates occurring at the time in the ab- 
scissas. We note that for each frequency the minimum contour level plotted 
is much greater than the signal uncertainty. We have overplotted the density 
flux at 15.4 GHz (dark thick curve), 212 GHz (light thick curve) and the time 
derivative of the GOES flux (black thin curve). The dashed curve represents 
the cone- of -influence of the transformation, i.e., where the boundaries do not 
affect the result. We observe that near peak time (11:44 UT) there are three 
groups of structures with scales around 10-20 seconds. The first one, a little 
before 11:44 UT, is observed at the three frequencies with no appreciable delay 
between them. The second one, just after 11:44 UT is present also at the three 
frequencies. The last one (~ 11:44:30 UT) is observed at 15.4 GHz and GOES 
time derivative. At 212 GHz, this structure appears only after 11:44:40 UT. 
Other structures with longer duration are also observed for the 40 seconds 
scale, around 11:45:50 UT and around 11:46:30 UT. In general, the wavelet 
transforms are similar for the different frequencies, including the coincidence 
in time of the characteristic features at the different scales. Time profiles do 
not reflect this in a more evident way because the relative importance of the 
pulses are different. 



3 The impulsive phase 



Figure 4 (left panel) depicts the time evolution of the 2.695 and 15.4 GHz 
emissions at the top, together with the temporal evolution of the optically 
thick spectral index at the bottom. This index is determined from the flux 
densities at 2.695, 4.995, 8.8 and 15.4 GHz. The index is lower than 1 all 
along the impulsive phase. The spectral index of the optically thick gyrosyn- 
chrotrqn spe ctrum originated in a homogeneous source is between 2.5 and 2.9 
( iDulkl . Il985l ). independently of the electron index. Therefore, our result is an 
indication of the presence of an inhomogeneous source. 



The 40 ms integrated time evolution of the flux density at 212 and 405 GHz 
during the impulsive phase is shown in the right panel of Figure 4 (above). The 
bottom of this panel corresponds to the temporal evolution of the optically 
thin spectral index deduced from these flux densities, the only two frequencies 
that are optically thin in our data. During the central part of the impulsive 
phase (between the dashed lines in the right panel of the Fig. 4), the index is 
close to constant, with gradual transitions at the beginning and end. Therefore, 
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Fig. 4. Left panel: 1 s time profiles of 2.695 and 15.4 GHz flux densities from RSTN 
observations (above) and temporal evolution of the optically thick spectral index 
a obtained from the 2.695, 4.995, 8.8 and 15.4 GHz flux densities (below). Right 
panel: 40 ms time profiles of 405 and 212 GHz flux densities deduced from SST 
observations (above) and temporal evolution of the optically thin spectral index a 
obtained from the 405 and 212 GHz flux densities ratios (below). 

we can characte rize the pea k time with only one electron population. Using 



Bulk's formula (IDulk 
(5^3. 



19851 ). we can estimate the electron energy index as 



The low spectral resolution of our radio data does not allow us to obtain a 
unique theoretical fit of the spectrum, i.e. different sets of parameters (mag- 
netic field strength, viewing angle, source size, lower and higher cutoff ener- 
gies, etc.) give similar results for the optically thin region; while, the optically 
thick region is never well fitted, as expected when an inhomogeneous source 
is present. 



4 The 212 GHz gradual phase 



The gradual phase is observed very clearly at 212 GHz (Fig. 5). At lower fre- 
quencies, after the main burst, new bursting episodes of gyroemission origin 
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Fig. 5. The 212 GHz gradual phase. After the burst (~ 11:45:30 UT), a gradual 
increase of the flux is observed at 212 GHz lasting until 11:51:30 UT. During the 
gradual phase, some bursting episodes are seen. In contrast, after 11:45:30 UT, the 
microwave flux, here represented by the shifted 15.4 GHz (dashed line), shows a 
second burst composed of many pulses. The most intense one is also observed at 
212 GHz. Along interval A we integrate the flux density at the different frequencies 
to obtain one spectrum (Fig. 6). The continuous thick line represents the expected 
emission of a thermal bremsstrahlung. 

are seen. Figure 5 shows the 15.4 GHz flux density, arbitrarily shifted by 300 
s.f.u., to clarify the plot; the late bursting episodes are evident in this figure. 
The strongest 15.4 GHz feature is also observed at 212 GHz (arrow in Fig. 5). 
Variations in the atmospheric attenuation and/or emission are relatively more 
intense at 405 GHz than at 212 GHz, and thus they modulate the background 
making the free-free emission difficult to observe. 

To test the thermal bremsstrahlung hypothesis we have modelled the observed 
emission as an isothermal homogeneous bremsstrahlung source. We integrate 
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Fig. 6. Spectrum of the flux density during the interval 11:48 - 11:49 UT. The solid 
curve represents the expected emission of an isothermal bremsstrahlung source. 
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Fig. 7. A soft X-ray image taken by SXI of the flaring area just before the maximum 
of the flare. The contours represent the emission levels of the thermal source. The 
circle corresponds to the source of our isothermal bremsstrahlung model. 
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the flux densities at different frequencies in the 1 minute interval shown in 
Figure 5 with a label 'A' to obtain the gradual spectrum. We have not tried 
to subtract the contribution from other emission mechanisms, hke gyrosyn- 
chrotron; therefore, we choose an interval that does not show strong gyroemis- 
sion. To determine the size of the sour c e we have used an i mage from the Solar 
X-ray Imager fGOES/SXl. iHill et all . [2OO5I : IPizzo et all l2005f ) (Fig. 7), con- 



sidering the emitting area above 30% of the maximum, which yields a source 
size of the order of 40" . The emission measure and temperature of the source 
are determined using the Chianti-based model with coronal abundances and 
the response of GOES detectors. We flt the expected emission to the observed 
data, having as free parameters the source size and the source depth and we 
include a multiplicative factor (. The time proflle of the gradual emission at 
212 GHz must be fltted simultaneously. We have obtained the best flt with 
a source size of 28", 8.4 10* cm depth and ( = 7.5. Figures 5 and 6 show 
the results. The temporal evolution of the model is remarkably similar to the 
gradual phase of the 212 GHz emission. Moreover, the modelled spectrum fol- 
lows the observed one. 



The homogeneous and isothermal thermal bremsstrahlung flux density, F^u) 
is: 



F(u) = C ^ ^""J (1 - exp {-K,H)) n 



where k^, is the absorption coefficient, H is the depth of the source and Q 
the source solid angle, while T represents the source temperature and ( is our 
multiphcatiy e factor. The absorption coefficient be approximated by 

(lDuMll985h : 



0.2 N?L J _i 



with Nmed the thermal plasma density. In the optically thick regime, where 
= H » 1, Eq. 1 tends to: 

- C 2 ^ • 3 



Therefore, if we want to explain the observed spectrum, the product TVL of 
the microwave source must be a factor 7.5 larger. In the optically thin regime, 
Ty « 1, Eq. 1 turns out to be: 
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0.4 kn (EM" 



(4) 



where EM = N^^^V is the emission measure, V the volume of the source 
and R is the Sun - Earth distance. From Eq. 4 we can infer that the prod- 
uct EM/\/T must be a factor 7.5 larger; i.e. either the EM is larger or the 
temperature is lower, or both things are true. Therefore, considering the opti- 
cally thick and optically thin regimes, we conclude that the source is neither 
homogeneous nor isoth ermal. Similar resu lts for the gradual phase of other 
ev ents were obtained bylLiithi et al.l (120041) at submillimeter wavelengths, and 
by IChertok et al.l (119951 ) and ISui et al.l (120051 ) in microwaves. 



5 Conclusions 



In this work we have presented a multiwavelength analysis of a GOES XI class 
flare, including radio data, SXR and EUV images. Radio observations cover 
the range between 1.4 and 405 GHz, but we have no observations between 15.4 
and 212 GHz. At 212 GHz, the event has two distinct phases: impulsive and 
gradual. At microwave frequencies, there is a second impulsive phase whose 
emission at submillimeter frequencies is barely observed. 



We find that the optically thin spectral index a is almost constant during the 
central part of the impulsive phase, afterwards there is a hardening. A transi- 
tion to free-free emission, that in the optically thin regime has an almost null 
index, can explain the change of a. 



The gradual phase observed at 212 GHz can be explained by bremsstrahlung 
from a thermal source. The similarity of the temporal evolution of the emission 
coming from an isothermal bremsstrahlung source (with EM and T deduced 
from SXR data) and the flux density at 212 GHz, sustains this hypothesis. 
Nonetheless, either the EM must be higher or T must be smaller to fit the 
observations. On the contrary, to fit to the microwave observations, the prod- 
uct TVt must be larger. Therefore, the thermal source cannot be homogeneous. 



The use of submillimeter frequencies in the analysis of burst phenomena opens 
new windows to understand the earlier stages of a flare at high energies, which 
are hard to observe with HXR detectors, and the late stages of the flare, where 
the thermal bremsstrahlung is normally optically thin and more intense than 
at microwave frequencies. 
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